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Abstract 

Three dimensional SO(3) gauged Skyrme models characterised by specific 
potentials imposing special asymptotic values on the chiral field are considered. 
These models are shown to support finite energy solutions with nonvanishing 
magnetic and electrix flux, whose energies are bounded from below by two 
distinct charges - the magnetic (monopole) charge and a non-integer version of 
the Baryon charge. Unit magnetic charge solutions are constructed numerically 
and their properties characterised by the chosen asymptotics and the Skyrme 
coupling are studied. For a particular value of the chosen asymptotics, charge- 
2 axially symmetric solutions are also constructed and the attractive nature of 
the like-monopoles of this system are exhibited. As an indication towards the 
possible existence of large clumps of monopoles, some consideration is given to 
axially symmetric monopoles of charges-2,3,4. 
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1 Introduction 



In the 0(4), or usual, Skyrme |T| model in 3 dimensions, the finite energy conditions 
specify the asymptotic value at large distances of the chiral field uniquely. Depending 
on the parametrisation, either the SU(2) valued field U or the S 3 valued field 0°, 
subject to |0 a | 2 , with a = 1,2,3,4, this asymptotic value is 

lim U = 1 or lim 4 = 1 . (1) 

r^oo r— >oo 

The fields U and </> a are related through 

U = <p a o a , = 0° a a , (2) 

where in terms of the Pauli matrices f, a a = (ir l , I) and a a = (—ir l , I). 

Often, the static Hamiltonian of the Skyrme system is augmented with a 'pion- 
mass' potential 

v;(0 a ) = 77w(i-0 4 ) , (3) 

consistent with the asymptotics ([[]). The only practical effect that the inclusion of this 
potential (|3|) has is, that it renders the asymptotic behaviour of the chiral function 
exponential, where in its absence this would have been a power decay. 

The situation is very different when the Skyrme model is gauged in one or 
other 0] gauging prescription, as a result of which the asymptotic value of the chiral 
field is not fixed uniquely by finite energy conditions. This feature of 3 dimensional 
gauged Skyrme models was considered and highlighted in ||. 

In the present work, we augment the 3 dimensional 50(3) gauged Skyrme model 
studied in || |J, with the potential^ 

V = A(coscj -0 4 ) 2 , n>iu>0, (4) 

whose effect is to specify the asymptotic values of the chiral field uniquely, consistent 
with finiteness of the energy. The new asymptotics are 

lim0 4 = — 1 lim 6 = cosu . (5) 

r^O r— >oo v ' 

It is clear from flj^) that the volume integral of the density that maps the field 
space to the configuration space, is not going to be an integer except in the case 
where u — 0. Thus the lower bounds labeled by this charge cannot be identified 
with the degree of the map, or the topological Baryon charge, except when uj = 0. 
Such a noniteger charge however does supply a legitimate lower bound on the energy 
integral. For want of a better name, we shall persist in calling such charges Qb(w), 
with the understanding that only Qb{0) is really the Baryon charge. 

In the generic case 7r > uo > 0, there will be an independent lower bound in 
addition to Qb(uj), namely the magnetic monopole flux which also depends 

1 Like (fj), this potential is also chosen such that it results in the exponential behaviour of (f> 4 
asymptotically. 
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on uj. These lower bounds will be stated explicitly below. The main feature of the 
dynamics characterised by the potential (f|), with uj ^ is that the SO (3) is broken 
down to U(l), with the residual Maxwell field described by the correponding ' Hooft- 
tensor supporting a magnetic flux. Like the charge Qb(uj), this magnetic flux fi(uj) 
is integer also only modulo a continuous factor depending on uj, and takes an integer 
value only when uj — f . The solutions we have found turn out, as expected, to respect 
both lower bounds n(ui) and Qb{oj). 

We have confirmed the existence of finite energy solutions bounded by the charges 
h(uj) and Qb{^) by numerical construction. An interesting result is that when uj ^ 0, 
the solution persists even when the Skyrme coupling constant vanishes, i.e. k 2 = in 
(0) and @ . This is not surprising since the presence of the Yang-Mills term satisfyies 
the (Derrick) scaling requirement independently of the Skyrme term, and in this case, 
the soliton is bounded from below only by the magnetic flux //. This will be explained 
in more detail in Section 2, where the model is defined and the two said lower bounds 
will be stated. Then in Section 3, we study the spherically symmetric solutions for 
various uj and k, and find the ranges of these two parameters for which solutions exist 
by numerical construction. In Section 4, we study the axially symmetric magnetic 
charge-2 solution for the particular value of the parameter uj — |, with a view to 
learning whether two like monopoles of that system can be in an attractive or a 
repulsive phase. We encounter the rather surprising result that even for n 2 = this 
is attractive, and then as expected it becomes even more attractive with increasing 
K2 > 0. Section 5 is devoted to summarising and discussing our results. 

2 The model and lower bounds 

The model is specified by the gauging prescription and is the 3 dimensional model 
used in ||, augmented by the potential (Q). Usually, we will treat this potential 
as a gedanken entity and will not exploit it save as an agency justifying the asymp- 
totics (§. In terms of the S 3 valued field a = (0 a ,0 4 ), a = 1,2,3, and the SO (3) 
gauge connection with curvature F", the covariant derivative is defined by the 
prescription 

D^ a = d^ a + e a ^Al<\>\ D^ = d^. (6) 

Since much of the analysis will be almost identical to that in (the relevant) Section 

3 of [Q, we will use the same notation here. This will enable us to present some of 
the new results without the necessity of repeating the detailed analyses leading to 
them. The model is described by the Lagrangian 

£ = -4\F« U \ 2 + ^?|£W>T - \4\D[,rD^ b \ 2 - V(0 4 ) (7) 
which in the temporal gauge Aq = yields the static Hamiltonian 

n = 4\F t f 2 + ^|A0T + \4Wd^ + v(<p A ) . (8) 
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The potential V((f) 4 ) in both (|7J) and (§) is that given by (||). In our study of the 
'monopole', we will be mainly concerned with the energy density functional (§), but 
we give the corresponding Lagrangian (0) too in anticipation of our discussion of the 
corresponding 'dyon' solution. 

We proceed to state the two distinct lower bounds on the energy, namely the vol- 
ume integral of (§). Both bounds, the 'magnetic monopole' charge and the noninteger 
'baryon charge', pertain to the generic asymptotics (§) with uj ^ 0, dictated by @. 

The first of these follows from the classic Bogomol'nyi inequality 

4\f*\ 2 + 7^iiA<n 2 > . (9) 

It is obvious that the left hand side of (|9]) can be replaced by Ti of (||), by adding 
suitable positive definite terms to it, resulting in 

H > KlK X e ijk d k {rF%) , (10) 

on the right hand side of which we recognise the U(l) 't Hooft-tensor, Q -Fj", of the 
residual gauge field responsible for the magnetic flux provided that uj ^ 0, so that 
|0 a | — > sinu 7^ asymptotically. 

To state the corresponding inequality for the other lower bound, we define the 
'baryon charge' density and its covariantised version, respectively, 

Qo = j^e^e^di^d^dk^ (11) 

Qg = ■^e ijk e abcd D i <l> a D j <l> b D k <l> c <j> d . (12) 

The volume integral of (|11]) is the noninteger 'baryon charge' 

d 3 xg = (tt - uj)N , (13) 



except when u = 0, when it is simply the integer N, the degree of the map or the 
usual baryon charge. The actual gauge invariant charge density which enters the 
relevant inequality is defined in terms of qg in (|12"D by p. [| 

Q = Qg + 3e ijk <p a F« D k ^ , (14) 

whose volume integral turns out to be equal to the 'baryon charge' (TIB]). 
As shown in ||, it follows that 

> ^ (15) 

" '1+9(^)4 

We can now conclude from (P) and (]T~5|) the two distinct lower bounds on the 
energy E = J d 3 x7i, namely the 'magnetic' and 'baryonic' lower bounds, following 
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from the asymptotics 



E > Akk^Ki smu (16) 

E > UnKlK > (tt — oS) . (17) 
- /l + 9 (^ 1 



The two inequalities (jig) and (|T7T) signal the possibility of finding finite energy 
solutions bounded from below, provided that the (Derrick) scaling requirement is 
satisfied, which for (§) in 3 dimensions, it is. For the limiting case of u = considered 
in ||, inequality (^) trivialises and ( JT7|) then coincides with the lower bound used in 
In the other limit when u = it, both (|T6|) and (|17D trivialise so we would expect 
to find no nontrivial solutions in this case. This will be confirmed by our numerical 
results to be given below. For generic values of uo between these two limiting values, 
both lower bounds are valid independently, and any nontrivial finite energy solution 
must respect these. This will also be confirmed by our numerical results. 

It should perhaps be pointed out that neither of the bounds (0) and (|l|) can be 
saturated. As we shall see in Section 3, for uj ^ finite energy solutions persist also 
for k,2 = 0. But even in that case, the inequality (|T0|) cannot be saturated. Thus in 
the K2 = model, we have a system which does not saturate a Bogomol'nyi bound 
and whose stress-energy tensor therefore never vanishes. It follows that the charge-2 
monopole of this model is either attractive or repulsive, a property which it shares 
with the usual (ungauged) Skyrme model [0], in the latter case as is well known it 
being attractive ||. We shall find in Section 4 that the model with tz 2 — supports 
solutions describing mutually attracting like monopoles. Then as expected, when the 
Skyrme coupling constant k 2 is switched on this binding energy will grow further as 
is usual with other theories involving Skyrme like kinetic terms. 

3 Spherically symmetric solutions 

This Section is divided into two Subsections, in the first of which we present the 
reduced one dimensional subsystems of (0) and (§), while in the second we present 
our numerical results. 

3.1 One dimensional subsystems 

As in 0, we impose the spherical symmetry thus 

Aq = K^g(r) x a , A?= , (18) 

(p a = sin/(r) x a , 4 = cos/(r). (19) 

The ensuing reduced one dimensional Lagrange and (static) Hamiltonian, subject 
to a suitable rescaling r — > x such that all constants but the Skyrme coupling 
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are suppressed, are respectively 



L — —2 (2a' 2 + ^U^) " \{* 2 f 2 + 2a 2 sin 2 /) - 2^ 2 sin 2 / (V + ^f^) 
+x 2 g' 2 + 2a 2 g 2 (20) 



H = 2 (2a' 2 + ^-^) +kx 2 f' 2 +2a 2 sin 2 /)+2 K a 2 sin 2 / ( 2f' 2 + ^^f] (2 l) 



Note that we have suppressed the potential terms arising from (^) in (p0|) and (|2lD , 
since as will be explained below, nearly all our numerical constructions will be carried 
out in the A = limit. 

While (|2l|) is positive definite, (p0| ) is not. The latter will be relevant only in some 
remarks below concerning the dyon solution of this system. 

Let us first consider the case of main interest, namely the monopole solutions 
of the equations following from the static energy density functional (0). The finite 
energy conditions require the following asymptotic values 

lim fix) = 71 , lim fix) = u (22) 

x^O V ' x^oo J v 1 v > 

lima(x) = 1 , lim a(x) = , (23) 

x^O v ' x^oo y ' v ' 

as long as u 7^ . This is the case of interest in the present work. (The to = 0, which 
is only a limiting case here, was studied in detail in 

The behaviours of the functions f(x) and a(x) in the iCl region are independent 
of the value of A in (f|) , and they are 

f(x) = 7i + F 1 x + o(x 3 ) , (x < 1) (24) 
a(x) = 1 + A x x 2 + o(x 4 ), (x<l) (25) 

In the x ^> 1 region, the asymptotic behaviour of the function a(x) is again 
independent of the value of A and is 

a(x) = A e -^ xsinuj , A = const. , (26) 

while that of the function f(x) does depend on A. In the limit of vanishing A and 
finite A, these are respectively the power and exponential decays 

/(*) = v + - + o(\) , (27) 

X \X Z J 



_ -V2Ax 

f{x) =00 + F , (28) 

x 

where F and F are constants which can be evaluated by the numerical process. 

Thus, like with the usual (ungauged) Skyrme model, the addition of this potential 
results in the exponential localisation of the chiral function f(r). In the numerical 
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work presented in the next two Subsections we have used A = throughout, since the 
qualitative properties of the solutions are unchanged when A > 0. This was verified 
in many typical cases and thereafter the potential (|) played a Gedanken role, rather 
like for the Prasad-Sommerfield limit of the Georgi-Glashow model, except that here 
the A = limit is not particularly interesting since it does not lead to the saturation 
of any Bogomol'nyi bound. 



3.1.1 Dyon solution 

Before proceeding to describe our numerical results concerning the monopole solutions 
of this model, we briefly allude to the corresponding dyon solutions. Following Julia 
and Zee || we vary the energy density (|(]) with respect to the functions a(r), g{r) 
and f{r). Since (^) is not positive definite, the radial Ansatz ( |l8|) and (|l9|) is not 
guaranteed to be consistent with the full Euler-Lagrange equations of this system. It 
has however been verified in || that Ql8l)-(|Ts|) is indeed consistent. 

The crucial equation that signals the existence of a dyon solution is the r ^> 1 
asymptotic equation arising from the variation of the function a(x) 

a" = a -g 2 + sin 2 / + (29) 

which in the x — > oo limit reduces to 

a" = (sin 2 uj — g 2 )a , 

which yields acceptable exponentially decaying solutions only when the asymptotic 
value of q = lim^oo g(x) satisfies the condition 

< q < sinu , (30) 

and otherwise leading to unacceptable oscillatory behaviour for the function a(r) 
asymptotically. 



With the asymptotic condition fl30|) , one has the behaviour 

g (x) = q h o(x~ 2 ) , (31) 

X 

in which the constant c is evaluated by the numerical integrations and parametrises 
the electric flux of the dyon. We do not repeat here the detailed results of of the 
numerical process as this is identical to that for the dyon of the Georgi-Glashow 
model, as presented in |§. The relevant analysis in Section V of Ref.|| can be 
adapted to the present model, by substituting the condition < q < 1 there, by (p0|). 

The only qualitative difference between the dyon of the Georgi-Glashow model 
and the dyon of the present model is, that unlike in the former case here there is no 
Prasad-Sommerfield limit. 
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3.2 Numerical results 



Solving the spherically symmetric equations for numerous values of the parameters 
k e]0,oo[ and uj e]0,7r[, strongly indicates that they admit at least one regular, 
finite-energy solution for each choice of these parameters. 

The behaviour of the solution in the limit uj —>■ is different according to the 
value of k and is strongly influenced by the pattern of solutions occuring in the case 
uj — 0. We briefly recall (see ||, [| for further details) that for to — 0, solutions with 
a(oo) = exist only for k > K cr , K cr m 0.697. For k g]0.0, 0.697[ the relevant solution 
has a(oo) = 1. 

Let us now describe how the solutions look like for fixed k, with uj varying. We 
have found it convenient to characterise the solutions by the value of the asymptotic 
coefficient F defined in Eq. (27). The evolution of this parameter is reported in Fig. 
1 for several values of k. For k > n cr the solutions are such that the function f(x) 
(resp. a(x)) decreases monotonically from 7r (resp. 1) for x = to uj (resp. 0) for 
x — ► oo. The parameter F is positive as seen in Fig. 1 for k = 1.0. The classical 
energy decreases monotonically when uj increases. In the limit uj — > the classical 
solution of || are smoothly approached. 

The behaviour of the solutions is more elaborate when k < K cr ; this is illustrated 
on Fig. 1 and on Fig. 2 for k — 0.1 . One new feature is that the parameter 
F undergoes a change of sign when uj varies from to it. For large enough values 
of uj (say, uj > uj max , uj max m 0.5 in the case k = 0.1), the profiles of f(x) and 
of a(x) monotonically decrease as functions of x and the classical energy decreases 
for uj increasing. For uj < uj max the function f(x) develops a local minimum at an 
intermediate value of x and the parameter F defined in Eq.(27) becomes negative. 
Moreover, the function a(x) develops a local minimum and a local maximum at finite 
values (say x m , xm) of the radial variable. In this region of uj the classical energy 
increases with uj. When the limit uj — > is considered our numerical analysis indicates 
that x m stays finite, xu increases and the value cl{xm) approaches a = 1 in such a 
way that the corresponding profile of the uj = solution is approached on [0,2m]. 
The corresponding value of the classical energy is also reproduced. These different 
features are illustrated on Figs. 1,2. 

Considered as a function of uj, we also observed (see Fig. 2) that the classical 
energy is maximal at an intermediate value of uj. Our numerical analysis indicates 
that the maximum is attained precisely for uj = uj max , i.e. at the value where the 
change of sign of the parameter F occurs. The ^-dependence of uj max is reported on 
Fig. 3. 

As further suggested by Fig. 1, in the region k ~ 0.5 , uj w 0.05 the pattern of 
the solutions becomes very complicated. We obtained strong numerical evidence that 
several branches of solutions exist in this region. That is to say e.g, that we find more 
than one solution for k = 0.5, uj = 0.05. However these new branches seem to exist 
on a very small domain of the parameter uj, the numerical analysis is therefore rather 
difficult in this region. Since the study of such details is not the aim of the present 
paper, we refrained from further pursuing our numerical analysis in this region. 
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To finish this section, we mention that, choosing k = 0, we were able to construct 
numerical solutions for uo > 1.1729. The analysis of the solutions in the limit k — > 
(with fixed uj < 1.1729), seems to lead to a discontinuity of the function f(x). 



4 Axially symmetric solutions 

This Section is divided in two Subsections as in the previous Section, in the first 
of which the axially symmetric Ansatz, and, the boundary conditions of the axially 
symmetric solutions, are stated. In the second Subsection, the numerical results are 
given. 

Our objective in this Section is to construct higher magnetic charge solutions, and 
in the first instance charge-2 axially symmetric solutions, with the aim of discovering 
whether like monopoles of this model are in an attractive or a repulsive phase. In 
this framework, we will restrict our analysis to monopole rather than dyon solutions, 
in the temporal gauge A = 0. 

The analysis carried out in this Section is less general than that given in the pre- 
vious Section for the spherically symmetric solutions. There, we studied the detailed 
dependence of the solutions on the parameter u specifying the dynamics. Having 
exposed these properties satisfactorily, we proceed to study the most natural subset 
of models here, namely those specified by u — |, supporting monopoles of integer 
magnetic charges. 

Within this u — | subset of models, we consider the models specified by the 
Skyrme coupling tt 2 , or, the effective parameter k for the range k > which includes 
interestingly the point k = 0. 



4.1 Ansatz and boundary conditions 

With magnetic charge, or azimuthal winding, n = 1,2,3,.., the axially symmetric 
Ansatz |1(J for the gauge field is 

1 r 



M 2r 



rl {Hxdr + (1 - H 2 ) rd6) - n (t?H 3 + r™ (1 - H A )) r sin 6dxf> , (32) 



and for the Skyrme field it is 

U = - (cos / II + sin / [(sin g sin 9 + cos g cos 6)t™ + (sin g cos 8 — cos g sin 0)tq\) . 

2 (33) 
where Hi, if 2, H3, H4, f and g are functions of the coordinates r and 9. The symbols 
r™, t'q and denote the dot products of the cartesian vector of Pauli matrices, 
t = (t x , T y , r z ), with the spatial unit vectors 



el 1 



(sin 9 cos ncf), sin 9 sin ncj), cos 9) 
(cos 9 cos ncf), cos 9 sin ncj), — sin 1 
(— sin ncf), cos ri(p, 0) , 



(34) 



9 



respectively. 

For n = 1, Hi = H 3 = 0, H 2 = = a(r), f = f(r) and g = 9 the spherically 
symmetric ansatz of (|18|)-(|i9|) is recovered. 

The residual U{1) gauge degree of freedom Jl0[ is fixed by the condition rd r H\ — 
dgH 2 = 0, which is just the Coulomb gauge in the two dimensional residual U{1) 
subsystem resulting from the imposition of radial symmetry in the x — y plane, i.e. 
with radius p = \fx 2 + y 2 . 

At the origin the boundary conditions for the gauge field functions read 

-f^2|r=0 = #4^=0 = 1; -f^l|r=0 = -f^3|r=0 = , (35) 

and for the Skyrme functions 

f\ r =o = 7T, d r g\ r=0 = . (36) 

For the gauge field to approach the asymptotic configuration of a monopole we 
choose 

-f^2|r=oo = -f^4|r=oo = 0, Hi\ r=OQ = Hs\ r=00 = , (37) 

and for the Skyrme field functions 

f\r=oo = g\ r =oo = & ■ (38) 

The boundary conditions along the p- and z-axis are determined by the symme- 
tries. For the gauge field functions symmetry considerations lead to the boundary 
conditions 

Hi\q =0 = -f^3|e=o = , dgH 2 \e=o = doH^\ e= Q = , 
Hi\e=w/2 = H 3 \ g=n / 2 = , dgH2\e =n /2 = de-f^i 1 0=71-/2 = 

along the well as the condition 

#2 1 0=0 = H± 1 0=0 . (40) 

Along these axes the Skyrme field functions satisfy the following boundary conditions; 

def\e=o = , deg\o=o = 1 , + \ 

9$f \e=n/2 = , d g g\o=n/2 = 1 



4.2 Numerical results 

Solving the set of partial differential equations numerically for the model characterised 
by w = 7r/2, and for different values of K, we find that even for k = there is 
only an attractive phase. As shown in Fig. 4, the difference 8E of the energy per 
winding number n between the n = 1 and n = 2 increases with increasing k. This 
is an indication that in this model 2- monopole bound states can exist. Although a 
definitive demonstration of this is well beyond the scope of the present work. (That 
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would involve finding the dependence of the interaction energy on the separation of 
the two monopoles.) 

Moreover, we calculated the energy E n (K), for different values of k and n. The 
values are given in the tables below. Unfortunately the numerical process becomes 
less reliable with increasing monopole charge n. As a result we restrict our numerical 
constructions to n < 4 only For k = 0, K = 3 and n = 5 we find 



n 


1 


2 


3 


4 


E n (0) 


2.95 


4.82 


6.63 


7.56 




n 


1 


2 


3 


4 


E n (3) 


3.40 


5.40 


7.30 


8.80 



n 


1 


2 


3 


4 


E n (5) 


3.48 


5.60 


7.50 


9.31 



^From this data we can deduce some quantitative information on the interaction 
energies of the monopoles leading to the possible formation of lumps of charges n < 4. 
To this end we define the following 'binding energy' corresponding to the energy 
needed to dissociate a charge-n lump into a charge-n — 1 and a charge-1 lump, divided 
by the energy of the charge-n lump. 



AEi n ~ 1,1 h 



[E w _i(k) + Ei(k)] -EM 
EJk) 



(42) 



The values for different k and n are given in the Table below. 



n 


A^"- 1 ' 1 >(0) 


A^ n - 1 ' 1 >(3) 


A^ n - 1 ' 1 >(5) 


2 


0.22 


0.26 


0.24 


3 


0.17 


0.21 


0.21 


4 


0.27 


0.22 


0.18 



This table shows that for all three values of k, AEj^ 1 ^' 1 ^ (k) remains positive, which 
is an indication for the existence of monopole lumps of charges up to n = 4. 

It is hard to extract any reliable conclusion from such meagre data, but the fact 
that the binding energies do not seem to decrease with increasing n is encouraging 
from the point of view of the possibility of finding very large monopole clumps. This 
question will be investigated elsewhere, using different numerical techniques. 



5 Summary and Discussion 

We have studied a particular variant of the Skyrme model gauged according to the 
prescription (||), equivalent to the commutator gauging with respect to the SU{2) 
gauge connection Ai = — ^Ai • r, 

D i U = d i U + [A i ,U] , (43) 
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which is augmented by the potential ([|). The function of this potential is to fix 
the boundary value of the Skyrme field at large distances which, unlike in the usual 
(ungauged) Skyrme model, is not fixed by the requirement of finite energy^. Thus we 
have considered a set of the gauged Skyrme models characterised by the parameter 
uj appearing in the potential (§|). In the usual Skyrme model uj = 0. 

The qualitative properties of the solitons, depending on the parameters uj and k 
characterising the models, are studied in the spherically symmetric case. 

The main effect of the boundary condition uj ^ is the breaking of the 5*0(3) 
gauge symmetry of the solution down to SO (2), asymptotically This is related to 
the nonvanishing VEV of the Skyrme field (fi a . This can be seen clearly from the 
second member of (|5]). In addition to this magnetic charge, a noninteger version 
of the Baryon number given by flTBp is also associated with this solution. We have 
verified the existence of finite energy solutions bounded from below by both these 
charges for the allowed ranges of uj. The latter (ranges) depend also on the (effective) 
Skyrme coupling k of the model. These ranges have been illustrated in Fig. 1. 

A surprising if not unexpected property of these models is, that the model char- 
acterised by k = does support a soliton. This could have been expected since in the 
presence of the Yang-Mills term, it is not necessary to have a Skyrme term to satisfy 
the (Derrick) scaling requirement. We found that the solitons of the n = models 
are generally quite similar to those of the the models with k^O, with one noticable 
qualitative difference. This concerns the restricted range of allowed uj for the k = 
model as seen in Fig. 1, (and in Fig. 3) which in addition to this qualitative feature 
also illustrates the fact that for small enough k the profile of the chiral function f(r) 
sinks below its asymptote uj and approaches it from below. 

By contrast when uj = 0, the model with k = cannot support a soliton because 
in that case the lower bound (|TDJ) disappears, leaving only the lower bound ( |T3"D in 
place, and the latter trivialises in the k = limit. This is also borne out by the 
graphs in Fig. 1. 

After exposing the qualitative features of the solitons in the spherically symmetric 
case, we studied axially symmetric solutions of the model. Our aim here was to 
discover if like monopoles of the uj ^ models are in attractive or repulsive phases. 
For this purpose we restricted ourselves to the u = | model which has the nice 
feature of having integer magnetic charge. We did however consider varying values 
of the (effective) Skyrme parameter n, including the distinguished case of k = 0. 

The main interest in the model specified by k = in this respect is that, like 
the usual (ungauged) Skyrme model, it exhibits no free coupling constants that can 
parametrise the crossover from an attractive to a repulsive phase. Like the latter it 
lacks also a neutral, or Bogomol'nyi saturated phase, where the (static Euclidean) 
stress-energy tensor would have vanished leading to noninteracting solitons. Thus 

2 Unlike here, in Refs.j^, an explicit potential was not employed. This is because there, only 
one value, u> = 0, was used so it was considered reasonable to treat the role of the corresponding 
potential in a putative capacity. This is justified since, as we have verified here too, the explicit 
use of a potential only affects the asyptotics of the chiral function /(r) and otherwise leads to no 
qualitative differences. 



12 



the unique phase in which the solitons are supported is of special interest, especially 
if it were attractive because then the switching on of k would most likely not result 
in a crossover to a repulsive phase. We have verified that this is precisely what 
happens, as illustrated in Fig. 4 for the magnetic charge-2 soliton. An outstanding 
problem in this context is the calculation of the (attractive) interaction energy of two 
1-monopoles as a function of their separation. This would demonstrate the existence 
of bound states definitively. 

In addition to verifying that these models support mutually attracting like 
monopoles, we sought some indications as to whether there is the possibility of forming 
bound states of monopole charge greater than 2. (Apart from its intrinsic interest, 
the formation of very large monopole clumps may be relevant in cosmology.) To this 
end, axially symmetric solutions with monopole charges n = 2,3,4, for the k = 0, 
k = 3, and the k = 5 models were constructed. It was seen that all these solitons 
remained in attractive phases. While it is expected that these axially symmetric 
solitons with n > 2 are not the lowest energy solutions, the latter probably exhibiting 



(solid) Platonic symmetries [|TT| [T2[ as in the usual Skyrme model, it is nonetheless 
true that they give a reliable indication towards the existence of such bound states. To 
this end we list the binding energies against the dissociation of an axially symmetric 
charge-n monopole into a charge-n — 1 and a charge- 1 monopole in the last Table 
of Section 4. We see that this binding energy stays positive and does not change too 
much as n increases, at least for n < 4. Unfortunately the numerical process was not 
reliable much beyond n = 4. 

While in the present work our primary aim has been the qualitative study of 
the solitons of the gauged Skyrme model charcterised by u> ^ 0, and the resulting 
property of the interaction of like monopoles in these theories, it may be in order to 
emphasise some physically attractive features of these models. These models combine 
some attractive features of (a) Higgs models, in this case the Georgi-Glashow model 
in that they support monopoles, and features of (b) Skyrme models, in this case the 
usual Skyrme model |]J in that the like-charged solitons are in an attractive phase. 
Indeed it appears that the attraction properties of these models are considerably more 
pronounced than those of the Skyrme model 

Apart from the practical consideration of the possibility of supporting large 
monopole clumps, there are two theoretical properties of the models that deserve 
mention. One is the physically desirable property of the symmetry breaking from 
SU{2) to U(l), and the other one is the fact that the k = model supports solitons 
with much the same qualitative properties as the generic models. The considerable 
advantage of this is that we avail of the Skyrme theortic feature of mutually attracting 
solitons without having to pay the price of featuring a (quartic kinetic) Skyrme term 
in the Lagrangian, thus avoiding the attendant severe problems of quantisation. 
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Figure 1: The quantity F, that determines the asymptotic behaviour of the skyrme 
field function / with f(x >> 1) = + F/x + o(l/x 2 ) is shown as a function of for 
n — 1 and different values of /c. 
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Figure 2: The energy E of the n = 1 solution is shown as a function of uj for « = 0.1 
and A = 0. The two inlets show the skyrme field function f(x) as function of the 
radial coordinate x for a) uj = 0.3 < uj max and b) uj — 0.9 > 0J max . 
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Figure 3: uj max , the value of uo at which the energy has its maximum, is shown as 
function of k. The inlet shows the quantitiy F (see Fig.l) over uj for k = and 0.1. 
The asteriks mark oj max {.^)- 
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Figure 4: The energy per winding number E/n is shown as a function of k for n — 1 
and n = 2,c<j = 7r/2 and A = 0. The inlet shows the difference of the energy per 
winding number SE = E(n = 1) — E(n = 2)/2 between the n = 1 and the n = 2 
solution. 
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